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^ . ABSTRACT 

The most massive elliptical galaxies apparently formed the fastest, because the ratio 
of a elements (such as oxygen) to iron is the smallest. In fact, iron is mainly pro- 
■ duced from type la supernovae on a timescale of ~ 0.1 — 1 billion years, while the a 

' elements come from massive stars on timescales of a few tens of million years (Mat- 

• , teucci 1994). Reproducing such a a/Fe correlation has long been a severe problem for 

i-^ ■ cosmological theories of galaxy formation, which envisage massive galaxies to assem- 

ble gradually from smaller progenitors, and to be characterized by a star formation 
Q . history too much extended towards late cosmic times. While it has recently become 

i-{ • clear that feedback from Active Galactic Nuclei (AGN) activity play a role in the late 

C/3 I quenching of star formation (e.g. Cattaneo et al. 2009), and that early star formation 

i history in the galaxy progenitors affect the a/Fe ratio (Calura & Menci 2009), major 

mergers alone cannot enhance the star formation in the high-redshift progenitors to 
the levels required to match the steepness of the observed a/Fe correlation (Spolaor 
^ ' et al. 2010). Here we report that the inclusion of the effects of fly- by 'harassments', 

00 ' that trigger lower level starbursts, combined with the AGN quenching of the starburst 

. activity, considerably enhances the capability to account for the observed a/Fe ratio 

0^ ' in ellipticals within cosmological galaxy formation models . The critical difference be- 

\l I tween the earlier work and the present result is the effect of starbursts driven by fly-by 

encounters that would have been very common amongst the high-redshift progenitors 
of massive galaxies and which would have boosted star formation in the first 2 billion 
years after the Big Bang, combined with quenching of the burst activity within the 
first 3-4 Gyr. 
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1 INTRODUCTION (Thomas et al. 2005; Nagashima et al. 2005; Spolaor et al. 

2010) in accounting for the above correlation is in fact due to 



The correlation between stellar a/Fe and velocity dispersion 
observed in local ellipticals (Trager et al. 2000; Thomas et 
al. 2005; Bernardi et al. 2006; Graves et al. 2007; Sanchez- 



such a "downsizing" aspect. In fact, the precision measure- 
ments of the power spectrum of the initial density perturba- 

, , „ , , „, , , tions (Percifal et al. 2007) provide evidence for large-scale 

Blazquez et al. 2006, Spolaor et al. 2010; Zhu et al. 2010) is n ^ \- ^ , „ „ 

, r, r, tiuctuations to be smaller, on average, compared to small- 

generally reierred to as one 01 the evidences 01 the downsiz- , . , , , , . 1 • ^ 1 1 1 

scale ones - m accord to what expected m a Cold Dark 

mg pattern or local galaxies, a pattern which is observable t,, .. /^tmt\ . . . 1. 1 . 1 . 1 1 .1 1 

Matter (CDM) scenario; this implies a hierarchical build- 
up to high redshift (e.g. Cowie et al. 2006; Kodama et al. r 1 1 /tmv/tn 11 1 n 1 ■ 

, , N , up ot dark matter (DM) haloes, where small-mass galaxies 

2004; Clements et al. 2008) indicating that the more massive , , 1, , r ■ n- .■ 1 

/ , collapse farst, and later assemble to torm massive ellipticals 

the galaxy, the shorter the star formation timescale (Mat- /t^ t ■ .1 mi , 1 • 1 ■ ■ r,i i- 

1 (lJ6 Lucia et al. 2006). the deep physical origin ot the dis- 
teucci 1994; Renzini 2006; Pipmo et al. 2010; Rogers et al. ^ 1 ^ 1 1 1 i- ^ 1 /r^ 

10 agreement between the observed and the predicted Q/i' e — (7 

correlation has been often considered as an evidence for a 
The lailure of cosmological galaxy formation models r ^ ^ 1 r 1 • 1 ^tmvj 11 r 1 

° o ^ fundamental flaw ot cosmological CDM models 01 galaxy 

formation (Thomas et al. 2005; Baugh 2006). 
* E-mail: fcalura@oats.inaf.it To produce downsizing within standard galaxy forma- 
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tion models, some physical mechanisms arc required which 
can enhance star formation in massive galaxies at early 
times, possibly quenching it at late times. In cosmological 
galaxy formation, environmental effects are likely to play a 
major role in producing such a scenario. The progenitors of 
massive galaxies formed from biased regions of the density 
field, and fly-by events and merging in such regions may trig- 
ger starbursts at early cosmic times {t < 2.5 Gyr), when the 
Universe was denser and the rate of galaxy encounters was 
at least 5-10 times higher than the present. A fundamen- 
tal mechanism to quench star formation in massive galaxies 
is linked to their Active Galactic Nuclei (AGN). Feedback 
from powerful AGN may switch off star formation when the 
progenitors assembled into the main progenitor of the final 
elliptical (Cattaneo et al. 2009). 

These processes may concur in yielding stellar popu- 
lations in massive galaxies characterized by a prompt star 
formation at early epochs followed by a quenching, in agree- 
ment with what suggested by the observed a/Fe—a correla- 
tion. Thus, testing whether the latter processes may in fact 
explain the above correlation is crucial for assessing the con- 
sistency of CDM models, and requires developed modelling 
which includes all the physical processed detailed above. 

In a previous paper (Calura & Menci 2009, hereinafter 
CM09), we investigated the chemical properties of local 
galaxies within a cosmological hierarchical clustering sce- 
nario through a semi-analytic model (SAM) of galaxy for- 
mation. We used a hierarchical semi-analytic model where 
chemical evolution was computed by taking into account the 
stellar lifetimes, a significant step forward with respect to 
the instantaneous recycling approximation. However, chem- 
ical abundances were computed by considering the total 
(summed over all progenitors) stax formation history of each 
model galaxy and assumed an average interstellar a/Fe. This 
treatment was suited to the study of interstellar abundances, 
comparable to the ones observed today in stars in the solar 
neighbourhood, but in many cases it lead us to underesti- 
mate the integrated stellar a/Fe abundances measured in 
local early type galaxies. 

In this paper, by means of the same SAM, we perform 
a thorough study of the the average stellar a/Fe in early 
type galaxies, taking into account the contribution from all 
the single progenitors. Although computationally more ex- 
pensive, this constitutes a step forward with respect to our 
earlier work (CM09), since thanks to the inclusion of the 
contribution of all the progenitors, our estimates of the in- 
tegrated abundances arc more accurate. A second enhance- 
ment with respect to our previous work is the calculation 
of the luminosity-weighted abundances (against the mass- 
weighted ones used in the first paper), another aspect which 
improves the comparison of our results to local observations. 
The aim of this paper is to show the relative roles of the var- 
ious processes in shaping the a/Fe — a correlation, stress- 
ing the importance of some processes which have been ne- 
glected in previous studies, such as interaction-triggered fly- 
by events and mcrging-triggered starbursts. 
This paper is organized as follows. In Sect. 2, our model 
is briefly described. In Sec. 3, we present our main results, 
which are then discussed in a broader context in Sect. 4. 



2 THE MODEL 

To study the integrated abundances of local ellipticals, we 
use the same model as described in Calura & Menci (2009). 
We defer the reader to that paper for further details. Here 
we summarize the main features of our model. 
Galaxy formation is driven by the collapse and growth of 
DM halos, which originate by gravitational instability of 
overdense regions in the primordial DM density field, a 
random, Gaussian density field with CDM power spectrum 
within the "concordance cosmology" (Spergel et al. 2007) 
for which we adopt round parameters Qa = 0.7, f2o = 0.3, 
baryonic density Q.^ — 0.04 and Hubble constant (in units 
of 100 km/s/Mpc) h = 0.7. As cosmic time increases, larger 
and larger regions of the density field collapse, and eventu- 
ally lead to the formation of groups and clusters of galaxies, 
which grow by merging with mass and redshift dependent 
rates provided by the Extended Press & Schechter formal- 
ism (Press & Schechter 1974). 

The clumps included into larger DM halos may survive as 
satellites, merge to form larger galaxies due to binary aggre- 
gations, or coalesce into the central dominant galaxy due to 
dynamical friction. The above complex dynamical evolution 
of galaxies is followed through a Monte Carlo simulation of 
the collapse and subsequent merging history of the peaks of 
the primordial density field, which enables to generate a syn- 
thetic catalogue of model galaxies and of their past merging 
history. 

For each model galaxy, we follow the evolution of its gas 
and stellar content; the gas which has radiatively cooled in 

the galactic haloes (with mass Mg) partly condenses into 
stars at a rate i/" oc Mg/td on long timescalcs, ~ 2 — 10 
Gyr according to the Kennicutt law (Kcnnicutt 1998) . Star 
formation can also occur through interaction-triggered star- 
bursts (ITS), driven by merging or by fly-by events between 
galaxies. Such a star formation mode can convert large frac- 
tions (up to 100 % in major merging events) of the cold gas 
on short timescales (ss 10® yrs), and provides an important 
contribution to the early formation of stars in massive galax- 
ies (Menci et al. 2004a). 

The model also includes a detailed treatment of the growth 
of supermassive black holes at the centre of galaxies, of 
the corresponding AGN activity powered by the interaction- 
triggered inflow of cold gas onto the black holes, and of the 
AGN feedback onto the galactic gas. Due to the crucial role 
of AGN feedback in the affecting the star formation history 
of massive galaxies, this sector of the model has been ex- 
tensive tested; the predicted QSO luminosity function and 
X-ray luminosity function and emissivity have been checked 
against observation over a wide range of redshift < z < 5 
and bolomctric luminosities 10*"^ < L/ergs~^ < 10*® (Menci 
et al. 2003; 2004a); the description of the AGN feedback as 
resulting from the sweeping of galactic gas by the blast wave 
produced by the injection of the AGN radiative energy has 
been tested by comparing the absorption properties of AGNs 
(as a function of both luminosity and redshift) with recent 
observation (Menci et al. 2008) ; the consistency of the mod- 
elling of the AGN feedback with the color distribution of 
galaxies has been discussed in Menci ct al. (2006). 
As shown in Fontana et al. (2004), the inclusion of pro- 
cesses such as ITS and AGN feedback within the SAM af- 
fects mostly the model rendition of the stellar mass function 
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at high rcdshift. All models discussed in this paper have ba- 
sically the same stellar mass function at z ~ 0. 
Galaixy luminosities in different bands are computed by con- 
volving the computed star formation history (SFIf) with 
synthetic spectral libraries of stellar populations (Bruzual 
& Chariot 2003). 

To compare our prediction to the observed a/Fe-a relation, 
as performed in our previous paper, we select elliptical galax- 
ies from the SAM synthetic catalogue on the basis of their 

prcscnt-day (B-V) colour. In specific, we consider elliptical 
galaxies all the systems with present (B-V)> 0.85 (Roberts 
& Haynes 1994). For any selected galaxy, we compute the 
chemical evolution a posteriori by means of detailed chem- 
ical evolution equations taking into account the stellar life- 
times, fundamental for studying the chemical evolution of 
Fe. Thus, the SFH, the gas content, and the gas inflow 
and outflow histories of all the progenitors of the consid- 
ered galaxy provided by the SAM constitute the input for 
calculating the chemical evolution. This constitutes a step 
forward with respect to our earlier work (CM09), where we 
considered only the total (summed over all progenitors) star 
formation history of each model galaxy and assumed an av- 
erage interstellar a/Fe. 



3 RESULTS 

The average stellar a/Fe may be calculated as: 

/ (Q/Fe)(t)7A(t)I.^^^(fo - t)dt 



< a/Fe >*: 



/ ^p{t)Lv(to - t)dt 



(1) 



i.e. it is the time integral of the interstellar (a/Fe) abun- 
dance over the star formation history, represented by the 
star formation rate '^{t), multiplied by the luminosity 
Ly^^{to —t), which is the luminosity at epoch to of a simple 
stellar population of unitary mass with an age of (to — t), 
and is given by SSP photometric models as a function of age 
and metallicity (BC03), divided by the present total lumi- 
nosity. 

In Fig. 1, we show the predicted a/Fe vs a relation for our 
selected elliptical galaxies, considering four different cases. 
In Fig. 1, we report two different sets of data (Arrigoni et 
al. 2009; Spolaor et al. 2010), in order to have an idea on 
the uncertainty of the slope of the observed a/Fe-a rela- 
tion. In the standard case, we consider a model including 
interaction-triggered starbursts in the most massive galaxies 
and Active Galactic Nuclei (AGN) feedback. In the second 
case, we suppressed the ITS and consider AGN feedback 
alone. In the third case, we suppressed AGN feedback and 
consider ITS alone. Since the ITS consist of two components 
(merging and fly-by events), it may be interesting to sepa^ 
rate these two processes in order to appreciate the impor- 
tance of fly-by events. In the fourth case, we show our results 
obtained having suppressed AGN feedback and fly-by inter- 
actions, but considering only merger-triggered starbursts. In 
the standard case, the slope of the observed a/Fe-a relation 
is satisfactorily reproduced. In the second and third cases, 
although with our model a correlation between a/Fe and a 
is obtained, the predicted relation is considerably shallower 
than the observed one. In the fourth case, the predicted 
a/Fe-cr relation is overly flat with respect to the observer 
tions. The comparison between the third and fourth cases is 



essential to understand the major role of fly-by interactions 
in driving the correlation between a/Fe and a. 



4 DISCUSSION 

Our results show that a standard CDM galaxy formation 

model including ITS, consisting in both merger-triggered 
and fly-by triggered starbursts, and AGN feedback can nat- 
urally lead to shorter star formation timescales in larger 
galaxies, i.e. to star formation histories reversed with re- 
spect to their mass-assembly histories. 

The absolute novel ingredient of our model consists in the 
simultaneous inclusion of both the starbursts triggered by 
fy-by events, absent in any previous investigation of this 
fundamental property of early-type galaxies, and the AGN 
feedback related to the active Quasar phase of AGNs. 
The general problem of quenching the star formation his- 
tories in massive galaxies have been investigated in several 
studies (e.g. Gabor et al. 2010) and is of great interest. Fig. 2 
is useful to better understand how the features of our model 
act directly on the star formation histories of single galax- 
ies. In each panel of this figure, we show the SFH of three 
galaxies drawn from the sample including ITS and AGN, 
including AGN but not ITS and from the one including ITS 
but not the effects of AGNs. The three galaxies shown in 
each panel present the same merger histories. It is worth to 
stress that in each case, the final stellar masses arc not the 
same, since the inclusion of physical ingredient such as ITS 
and AGN can considerably alter the star formation hisory 
and the feedback history of the systems. The final stellar 
mass of each model is reported in Fig. 2 (see caption for 
further details). 

It is important to note that the suppression of ITS has effects 
also on the efficiency of the AGN feedback, since, according 
to the model used here and described in Mcnci et al. (2008), 
the cold gas which causes the interact-triggered starbursts 
feeds also the nuclear activity. In fact, in every panel, both 
models without ITS and AGN present higher initial star 
fromation rates with respect to the standard case. The ex- 
planation of this is that in both cases (no ITS and no AGN) 
the AGN is inefficient in heating and ejecting gas from the 
galaxies, resulting in higher amount of cold gas when the 
star formation has already turned on. 

At the present time, the models which include ITS but with- 
out AGN feedback show SFR values larger than those shown 
by the ones drawn from the standard sample. Globally, their 
SFHs are less skewed towards early times with respect to 
the ones of the model with both AGN and ITS. The effects 
of AGNs in quenching the star formation activity within 
the first 3-4 Gyr is evident from this figure. Finally, the 
model with AGN but without ITS presents very extended 
SFHs, with rather flat and smooth behaviours and very high 
present-day values. Without the effect of ITS, gas consump- 
tion at early times is not very efficient and at later times, 
after a few Gyrs, the large amounts of cold star-forming gas 
prevent AGNs from being effective and quenching the star 
formation. It is therefore clear that, in order to have realistic 
star formation histories for early type galaxies, both effects 
of ITS and AGN feedback are necessary. 
Both effects of AGN feedback and ITS seem to play an im- 
portant role in determining galactic downsizing. However, 
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the effect of the former seems shghtly dominant, and this 
is supported by the fact that the predicted a/Fe-a relation 
without ITS is sUghtly shallower than the one without AGN. 
Our study of the a/Fe-a relation tells us that our galaxy for- 
mation model can naturally account for one of the most im- 
portant aspects of the "downsizing" character of early type 
galaxies. Downsizing (i.e. oldest stars are in most massive 
galaxies, star formation shifts to lower masses at late times) 
in hierarchical models will arise if there is a minimum mass 
required to have star formation and a maximum mass above 
which feedback from AGN suppresses further star formation 
(Sheth 2003; Croton & Farrar 2008). Together, these will 
also lead to an [a/FcJ-sigma relation, since in principle, in 
such a scenario, the distribution of the formation timcscales 
will be narrower for the stars forming in the most massive 
systems. In this framework, ITS help because it allows to 
use up more of the gas earlier, thus leaving less work to 
be done by the AGN feedback, and leading to a steeper 
[alpha/Fe]-sigma relation. A detailed calculation of chemi- 
cal abundances in more hirarchical models (such as the ones 
above) will be of further help in understanding the relevance 
of these processes with respect to the downsizing character 
of local galaxies. 

It is important to stress that the observations can be repro- 
duced without the modification of any fundamental chemi- 
cal evolution parameter, such as the stellar IMF. However, 
it is worth to stress that such ingredient may play some 
role in dctcrming the shape of the a/Fe-a relation. In fact. 
Although our unprecedented rendition of this relation, the 
predicted slope obtained with our standard model is slightly 
shallower than the observed one. This aspect will be inves- 
tigated in our future work, and may have its explanation 
in effects such as a possible IMF dependence on the SFR 
and/or a flatter IMF in starbursts (Recchi et al. 2009; Haas 
& Anders 2010; Calura et al. 2010). Some direct evidences in 
favour of the latter hypotcsis come from observational stud- 
ies of nuclear star clusters in the Milky Way (Bartko et al. 
2010). Also dynamical investigations of early- type galaxies 
seem to require stronger baryonic feedback at the epoch of 
the formation of these systems, fully consistent with an IMF 
skewed towards more massive stars (Napolitano et al. 2010), 
or from the study of local ultra-compact galaxies (Dabring- 
hausen et al. 2009). A top-heavy IMF at early times seems 
also required in order to reproduce the abundance ratios ob- 
served in the hot intracluster medium (Matteucci & Gibson 
1995; Locwcnstein & Mushotzky 1996). 

In our previous work (C09), the average stellar abundance 
ratios have been calculated without taking into account the 
contribution of the single progenitors, i.e. by considering the 
total star formation history of the selected galaxy, given by 
the sum of the star formation rates of all the progenitors, 
and considering an average interstellar a/Fe. In that paper, 
the observed a-Fe-cr relation could be partially explained 
by assuming a top-heavy IMF in the most massive galaxies, 
which mimics the effect of a higher star-formation efficiency 
(or, in other words, a shorter star formation timcscalc) in 
largest galaxies (Matteucci 1994; Ferreras & Silk 2003; Mat- 
teucci et al. 1998). 

The largest galaxies have the most massive progenitors, 
which form all their stars in the shortest timcscales and dom- 
inate the integral of Eq. 1. This is the main reason why here, 
by considering the SFRs and the abundance ratios of each 



single progenitor, instead of a total SFR and an abundance 
averaged over all the progenitors, we can successfully repro- 
duce the observed a/Fe-a relation without any modification 
of the main chemical evolution parameters. 
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Figure 1. Stellar average a/Fe ratio vs velocity dispersion compared to local observations. The colour code represents the predicted 
number of galaxies with a given a/Fe and a velocity dispersion a, normalised to the total number of galaxies with that velocity dispersion. 
The solid lines represent linear-regression fits to the model. The dashed line and dash-dotted line represent linear-regression fits to the 
observational a/Fe-a relations, as compiled by Spolaor et al. (2010; open diamonds) and Arrigoni et al. (2009; solid circles) respectively. 
Top-left panel: our results have been computed by means of our standard assumptions, i.e. by taking into account both interaction- 
triggered starbusts at high redshift and AGN feedback. Top-right panel: no interaction-triggered starbursts and fiy-by interactions. 
Bottom-left panel: no AGN feedback. Bottom-right panel: no AGN feedback, no fiy-by interactions. 
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Figure 2. Star Formation histories of a few early-type galaxies of our model. In each panel, we show the SFH of three galaxies drawn 
from the sample including ITS and AGN (solid black lines), including AGN but not ITS (cyan dotted lines) and including ITS but not 
AGNs (dashed red slines). In each panel, the three systems present the same merging history. The final stellar masses of the models are 
indicated in each panel. From top to bottom, the indicated masses refer to the standard case, the No ITS case and the No AGN case. 
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